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Abstract 



We analyze the effect of the Hall term in the magnetohydrodynamic turbulence under a strong 
externally supported magnetic field, seeing how this changes the energy cascade, the characteristic 
:ales of the flow and the dynamics of global magnitudes, with particular interest in the dissipation. 
Numerical simulations of freely evolving three-dimensional reduced magnetohydrodynamics 



q | scales of the flow and the dynamics of global magnitudes, with particular interest in the dissipation. 

•i-H 

O ■ (RHMHD) are performed, for different values of the Hall parameter (the ratio of the ion skin 



depth to the macroscopic scale of the turbulence) controlling the impact of the Hall term. The 
Hall effect modifies the transfer of energy across scales, slowing down the transfer of energy from 
the large scales up to the Hall scale (ion skin depth) and carrying faster the energy from the Hall 
scale to smaller scales. The final outcome is an effective shift of the dissipation scale to larger 
scales but also a development of smaller scales. Current sheets (fundamental structures for energy 
dissipation) are affected in two ways by increasing the Hall effect, with a widening but at the same 
time generating an internal structure within them. In the case where the Hall term is sufficiently 



X 

C3 , intense, the current sheet is fully delocalized. The effect appears to reduce impulsive effects in the 

flow, making it less intermittent. 

PACS numbers: 



I. INTRODUCTION 



Amone, various kinetic corrections to magnetohydrodinamics models (MHD), the Hall 
effect |l|, |2| has been considered of particular importance in numerous studies: magnetic 



reconnection 



regimes 



11 



lQJ? 



3], dynamo mechanisms 



accretion disks 



mm 



and physics of turbulent 



14J are some of the main examples. In this paper we want to study the general 



effect of the Hall term in magnetohydrodynamic turbulence in plasmas embedded in a strong 
uniform magnetic field, through numerical simulations. We studied the effect of this term 
on the dynamics of global magnitudes, the cascade of energy, the characteristic scales and 
the intermittency of the flow. 

The MHD models (one-fluid models) are important frameworks for the understanding 
of the large scale dynamics of a plasma. However, these models fail to describe plasma 
phenomena with characteristic length scales smaller than the ion skin depth pa = c/u P i 
(with ujpi the ion plasma frequency and c speed of light). At this level, the Hall effect, 
which takes into account the separation between electrons and ions, becomes relevant. To 
describe this regime is common to use the Hall MHD approximation, wich considers two- 
fluid effects through a generalized Ohm's law which includes the Hall current. In presence 
of a strong external magnetic field a new reduced model has been proposed, the RHMHD 
model 15Ml7J|). In this approximation, the fast compressional Alfven mode is eliminated, 
while the shear Alfven and the slow magnetosonic modes are retained |l8(. This new model 
(RHMHD) is an extension (including the Hall effect) of the previously known reduced MHD 
(RMHD) model. The RMHD equations have been used to investigate a variety of problems 
such as current sheet formation |l9|, |20(, non-stationary reconnection 2l|, |22j, the dynamics 



of coronal loops 23, |24|, and the development of turbulence 25j. The self-consistency of 
the RMHD approximation has been analyzed in ref. 26[. Moreover, numerical simulations 
have studied the validity of the RMHD equations by directly comparing its predictions with 
the compressible MHD equations in a turbulent regime 23]. In the same way it has been 
studied the validity of the RHMHD model (l5|. 

The control parameter (the Hall parameter) in this regime is e = pa/L, the ratio of the 
ion skin depth pa to the characteristic (large) scale of the turbulence L. The influence of 
the Hall term can then be studied by increasing the Hall parameter in different simulations. 

The organization of the paper is as follows: Section [Til describes the sets of equations 



used, and the codes to numerically integrate these equations. In Section |TTT] we present the 
numerical results: first a comparison between simulations with different Hall parameter is 
performed for some global magnitudes, second we study the energy spectra, and then we use 
different techniques to see the evolution of the characteristic structures of the flow. Finally, 
in Section IIVI we list our conclusions. 

II. EQUATIONS AND NUMERICAL SIMULATIONS 

The compressible Hall MHD equations (dimensionless version) are 



dV Tr 1 J x B „ (V 2 p^ 1 

V X W + — r V ' 



dt M\ p \2 M|( 7 - 1) 

V 2 V 1 V(V-V) 



.- p +(c + r )^ 7 -^, a) 
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V ■ A = (4) 

where V is the velocity field, w is the vorticity, J is the current, B the magnetic field, p 
is the density of plasma, A and <fi are the magnetic and electric potential. As indicated in 
Eq. ([1]) we used the barotropic law for the fluid, p = cte.p 1 (here p is the pressure), in our 
case we consider 7 = 5/3. Ms is the sonic Mach number, Ma is the Alfven Mach number, 
v and ( are the viscosities , 77 is the resistivity and e = pu/L (pa is the ion skin depth 
and L the characteristic scale of the turbulence) the Hall coefficient. All these numbers are 
control parameters in the numerical simulations. The Hall parameter e appears in front of 
the Hall term in the dimensionless equations, expressing the fact that the Hall term becomes 
important at scales smaller than the ion skin depth pa. 



A. RHMHD model 



The RHMHD model derived in the work by Gomez et al 16[ and numerically tested 



by Martin et al (15| is a description of the two-fluid plasma dynamics in a strong external 
magnetic field. The model assumes that the normalized (and dimensionless) magnetic field 
is of the form (the external field is along e z ) 

B = S Z + SB, |£B|«a<l, (5) 

where a represents the typical tilt of magnetic field lines with respect to the e z direction, 
thus one expects 

Vi«l, d x « a < 1. (6) 

To ensure that the magnetic field B remains divergence free, it is assumed that 

B = e z + \7 x(ae z + ge x ). (7) 

The velocity field, in the more general case, can be decomposed as a superposition of a 
solenoidal part (incompressible flow) plus the gradient of a scalar field (irrotational flow), 
i.e. 

V = V x (<pe z + fe x ) + V^, (8) 

where the potentials a(r,t), g(r,t), ip(r,t) and f{r,t) are all assumed of order a C 1 and 



ijj(r,t) is of order a expressing the possibility a slight compressibility (see details in [15 



17J). Introducing the expressions ([7j) and ([8]) in the compressible set of equations (|Tj)-([4j) 



du dj 



and taking the terms up to first and second order in a the RHMHD model is obtained: 

, [j,a]-[u;,^ + z/V 2 u;, (9) 

[ip, a] - e[b, a] + rjV 2 a, (10) 



da d(cp — eb) 
dt dz 



db n d(u - ej) r 7l n . . 
dt = P p dz + ^' ^ + ^ p[n ' ^ + 



-e^ p [j,a]+f3 p r ] V 2 b, (11) 



— = — + [<p,u]-[a,b]+vV 2 u, (12) 



where 



w = -V 2 ± <p, (13) 

J = -Via, (14) 

b = -d y g, (15) 

u = -d y f, (16) 

and the notation [A, B] = d x Ad y B — d x Bd y A is employed. (3 P = (3j/(l + (3j) is a function 
of the plasma "feeta" . 

We use this set of equations to study how the Hall effect modifies the dynamics of mag- 
netohydrodynamic turbulence under a strong magnetic field. 

B. Numerical codes 

We use a pseudospectral code to solve the set of equations (I9l- (fl~2l . Periodic boundary 
conditions are assumed in all directions of a cube of side 2iiL (where L ~ 1 is the initial 
correlation length of the fluctuations, defined as the length unit). In the codes, Fourier 
components of the fluctuations are evolved in time, starting from a specified set of Fourier 
modes (see section UTTl for the specific initial conditions), with given total energy and random 
phases. 

The same resolution is used in all simulations, 512 2 in the perpendicular directions to 
the external magnetic field and 32 in the parallel direction (this is possible because the 
structures that require high resolution only take place in the directions perpendicular to 
the field), allowing four different runs to be done with four different Hall coeficients. The 
kinetic and magnetic Reynolds numbers are defined as R = 1/V, R m = 1/rj, based on unit 
initial rms velocity fluctuation, unit length and non-dimensional values for the viscosity and 
diffusivity. Here we used R = R m = 1600 (u — ( — 1/1600, 77 = 1/1600) in all runs. We 
also considered a Mach number Ms = 1/4, Alfven number Ma = 1, and the Hall coeficients 
e=0, 1/32, 1/16, and 1/8 in all runs. 



A second-order Runge-Kutta time integration is performed, the nonlinear terms are eval 



uated using the standard pseudospectral procedure |28|. The runs are freely evolved for 
10 time units (the initial eddy turnover time is defined in terms of the initial rms velocity 
fluctuation and unit length). The magnetic field fluctuations were less than ten percent of 
the external magnetic field value, so we are in the range of validity of the RHMHD model. 

III. RESULTS AND DISCUSSION 

We performed simulations for four different Hall coeficients, e=0, 1/32, 1/16, and 1/8. 

To generate the initial conditions, we consider initial Fourier modes (for magnetic and 
velocity field fluctuations) in a shell in k-space 1 < k < 2 at low wavenumbers, with constant 
amplitude and random phases. Only plane-polarized fluctuations (transverse to the mean 
magnetic field) are included, so these are (low- to high-frequency) Alfven mode fluctuations 
and not magnetosonic modes. 

The runs performed throughout this paper do not contain any magnetic or velocity stirring 
terms, so the RHMHD system evolves freely. 

We study the influence of the Hall term in global quantities associated with the dissi- 
pation. Figures [1] and |2] show the mean square current density < J 2 > and mean square 
vorticity < u 2 > as function of time for e=0, 1/32, 1/16, 1/8. Both < J 2 > and < uj 2 > 
show that as the Hall parameter is increased the dissipation decreases (in the case of mean 
square vorticity this effect is considerably larger). Another remarkable effect is the shift in 
the peaks of these functions: < J 2 > and < u 2 > take longer to reach its maximum with 
increasing e. The time of the peak indicates the time where all spatial scales were developed 
(and therefore turbulence is fully developed). 

Figure [3] shows < J 2 + u> 2 > as function of time, the difference between the peaks is more 
clear in this case. Here we see two effects that occur simultaneously as the Hall coefficient 
is increased: The decrease in the dissipation and the delay in reaching the maximum point 
(and hence the time that it takes to develop all the scales). The first effect will have a direct 
impact on the dissipation scale of the respective flows while the second shows how the Hall 
term modifies its characteristic times. 

It is relevant to note that the dissipation scale (l/K^iss) is related to the number of scales 
that develop in the flow. It is common to consider that the decrease in the dissipation scale 
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FIG. 1: Current density, < J 2 >, as function of time for e = 0, 1/32, 1/16 and 1/8. The color of 
these curves are, black, violet, red and green respectively (color online). The thickness of the lines 
decreases proportionally to the value of e, thus, the thicker line corresponds to e = and the finest 
line to e = 1/8. The vertical straight line indicates a particular time where all the scales have been 
developed in all runs. Besides, in this time the value of < J 2 > is approximately the same for all 
runs. This particular time will be used to study the different structures in the flows. 



increases the range of developed scales in the flow (usually increases the size of the inertial 
range). In the same manner, the increase in the scale of dissipation leads to a decrease in 
the number of scales developed in the flow. However, it is not always this case. The results 
that we will show below indicate that the Hall term affects the total width of the dissipation 
range decreasing mildly the Kd iss (and therefore mildly increasing the dissipation scale) with 
the increase of the e, at the same time the delays suffered by the dissipation peaks is due 
to the development of a greater number of scales in the dissipative range due to a major 
accumulation of energy in these scales. 

To quantify the dissipation scale (in Fourier space) of the different flows we use the 



conventional criteria 



29| given by the equation ( ITTjl . 



K, 



< u 2 > + < J 2 > 



2 ^\l/4 



(17) 



In Table [J the Hall scale is shown along with the dissipation scale for each one of the 




15 



time 



FIG. 2: Vorticity, < uj 2 >, as function of time for e = 0, 1/32, 1/16 and 1/8. We are using the 
same convention of color /thickness than in figure 1 
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FIG. 3: < J 2 + uj 2 > as function of time for e = 0, 1/32, 1/16 and 1/8. We are using the same 
convention of color /thickness than in figure 1. The vertical straight lines indicate the maximum 
value for each curve. 





TABLE I: Hall and dissipation scales for the different runs. 




Run 


e KhclU 


-ft- diss 



1 132.25 

2 1/32 32 125.40 

3 1/16 16 124.58 

4 1/8 8 120.08 

flows. Here we see the decrease of the Kdiss in quantitative form with the increase of the 
Hall coefficient. Note that K m.** < K max = N/3 = 170 means that the runs are marginally 



31| for more demanding requirements if higher order statistic 



resolved (see Wan et al 2010 
is performed) 

The decrease of the global dissipation with the Hall parameter and the increase in the time 
of the peak development can not be understood by looking only at the temporal evolution 
of the global magnitudes. These effects could be due to a change in the characteristic time 
of the energy flow or to the development of small scale structures. 

To better understand these issues we study the energy spectra and the size and shape of 
the structures generated in the four runs. This help us to see whether or not the Hall effect 
produces the development of small scales and also to understand this dynamic in terms of 
the energy distribution. 

Looking at the spectra we can see the distribution of the energy through different scales. 
Figure |4] compares the energy spectra for all runs and Figure [5] shows a zoom around of the 
Hall scales used. 

As the Hall parameter is increased the energy spectrum is steeper at intermediate scales 
preceeding the dissipation range. At the same time there is an increase in the energy on 
scales smaller (larger k) than the dissipation scale (see Figs. |4]and[5]). The effect of the Hall 
term is then twofold: first there is a slow down of the energy transfer up to the Hall scale, 
resulting in a steeper spectrum, and then there seems to be a driving of energy from the Hall 
scale up to the small scales (see [30| for a study of how the Hall term affects the transfer of 
energy at different scales). A shift of the effective dissipation scale to larger scales is then 
to be expected (as indicated by the values of K,n ss given before) as well as a decrease in 
the global dissipation values. At the same time, since the Hall term increases the number 
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FIG. 4: Energy spectra for e = (solid line), e = 1/32 (dash-dotted), e = 1/16 (dash) and 
e = 1/8 (dotted). The vertical straight lines indicate the differents values of Kjjaii = V e f° r 
e = 1/32, 1/16, 1/8. The vertical straight dashed lines show the minimun and the maximum values 
of the Kdiss ( e = 1/8 and e = 0), through them we can see the effect of the Hall term on the 
dissipation scale. 

of effective scales on which the dynamics occurs (as evindenced by the extended spectra at 
small scales) a longer time to reach the peak of dissipation is expected, as previously shown. 

We study the characteristic structures of the flow and the effect of the Hall term by looking 
at the current density field. Figures [Ml] show the parallel component of the current density 
in a perpendicular plane to the external magnetic field at a given time for the different runs. 
The time was chosen in which all scales have been developed for all the flows (this time is 
indicated in Figured]). Also, for this particular time, the value of < J 2 > is approximately 
the same for all the runs. 

In Figure [6] (e = 0) we can clearly distinguish the current sheets that form in the flow. We 
have highlighted one of the current sheets with a rectangle with dashed lines. This structure 
is localized and well defined. Looking at the change of this structure with the value of Hall 
parameter we can see two effects: first, a widening of the sheet and secondly an internal 
filamentation. The widening is very clear from Figure O with e = to Figure [7] with e = 1/32 
and the internal filamentation starts to be seen in the Figure [BJ with e = 1/16, where also 
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FIG. 5: Enlarged view of a section of the energy spectra (limited k). The vertical straight lines 
indicate the different values of Knaii = l/ e for e = 1/32, 1/16, 1/8, these lines intersect each of the 
corresponding curves. 




FIG. 6: Parallel component of the current density in a perpendicular plane to the external magnetic 
field in the case with e = .Tones indicate out of plane current, with light tones=positive and dark 
tones=negative. 
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FIG. 7: Parallel component of the current density in a perpendicular plane to the external magnetic 
field in the case with e = 1/32 .Tones indicate out of plane current, with light tones=positive and 
dark tones=negative. 




FIG. 8: Parallel component of the current density in a perpendicular plane to the external magnetic 
field in the case with e = 1/16 .Tones indicate out of plane current, with light tones=positive and 
dark tones=negative. 
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FIG. 9: Parallel component of the current density in a perpendicular plane to the external magnetic 
field in the case with e = 1/8 .Tones indicate out of plane current, with light tones=positive and 
dark tones=negative. 

the thickness has increased. In the case with higher e = 1/8 the current sheet is completely 
filamentated, and is hard to distinguish a clear structure at all. 

These results are complementary to the results observed in the spectra and global mag- 
nitudes and corroborate the idea that the Hall effect results in an effective shift of the dis- 
sipation scale (current sheet thickness getting larger) but also an increase in the dynamical 
scale range (increase of filamentation). 

To better quantify the effect we have just observed, we plot the profile of the current 
density in the direction perpendicular to the current sheet seen in the figures EH These 
profiles are shown in Figure ITUl 

The net flow of current (the absolute value) is the same within the clear lines (vertical 
outside lines). When e = the current sheet is perfectly located (the dark lines mark the 
original position of the current sheet when e = 0) and it is homogeneous (in the sense that 
we have a single well defined peak). When e = 1/32 the original sheet expands and two 
sheets or filaments appear in their place (there are now two peaks). For e = 1/16 the width 
of the main sheet is greater and there is now a clear internal structure. In this case the 
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FIG. 10: Current density profile of the current sheets studied. From top to bottom are the cases 
with e = 0, e = 1/32, e = 1/16, and e = 1/8. The net flow of current (the absolute value) is the 
same within the clear lines (vertical outside lines), and the dark lines mark the original position of 
the current sheet when e = 0. 

ambiguity that arises is whether we have one or more sheets of current (compare Figure HD1 
with [8]) and hence the ambiguity of whether we have a wider sheet or two thin sheets. When 
e = 1/8 there is no trace of the current sheet. 

At this point we should make an important observation about the evolution of current 
sheets as a function of the Hall parameter. As we saw there are two effects acting simulta- 
neously, the widening of what could be considered the overall structure of the sheet and the 
internal filamentation that this suffers. In this way it could be interpreted that the Hall effect 
widens the current sheets (if we see the entire structure like the sheet) or on the other hand 
the Hall effect produces finer sheets (considering that the small filaments are the sheets). To 
remove the ambiguity (in semantics), we propose to speak in terms of dissipation, so if the 
global structure dissipates less energy as we increase the Hall parameter we will say that the 
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sheet is being widened, otherwise, if more energy is dissipated we will say that the relation 
between size and intensity of internal filaments allow us to identify new current sheets. Our 
results agree with the first frame of mind: as a function of dissipation the current sheets 
are widening and even more when e = 1/8 there is no trace of any structure that could be 
identified as a current sheet. 

IV. CONCLUSIONS 

We performed numerical simulations of magnetohydrodynamic turbulence in strong mag- 
netic fields, including the Hall effect, and varying the Hall parameter. 

We found that the Hall term affects the scales that are situated between the Hall scale 
and the dissipation scale, resulting in a decrease in the accumulation of energy in this scale 
range. The result is an effective shift of the dissipation scale but also a transfer of energy 
to smaller scales. When the separation between the Hall scale and the dissipation scale is 
larger an increasingly sharp steepening of the energy spectrum occurs at this range of scales. 
The final outcome is the generation of smaller scales when the Hall scale increases. 

Localized structures are destroyed by this effect, suffering a gradual filamentation with 
the increase of the Hall scale. The latter effect is manifested, for example, in the widening 
of the current sheets and the formation of internal structures within the sheets. At the same 
time a decrease of the total energy dissipated is observed. 

The results presented here suggest that the Hall effect reduces the intermittency, however 
a more detailed study of this property should be performed. We defer this to further work. 
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